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It is proposed that the state space of a quantum object with a complicated discrete 
spectrum can be used as a basis for multiqubit recording and processing of information in a 
quantum computer. As an example, nuclear spin 3/2 is considered. The possibilities of 
writing and reading two quantum bits of information, preparation of the initial state, 
implementation of the "rotation" and "controlled negation" operations, which are sufficient 
for constructing any algorithms, are demonstrated. 

We shall consider a nucleus with spin I = 3/2 and with an electric quadrupole 
moment in an axial symmetric crystalline electric field and a constant magnetic field 
parallel to the symmetry axis. We assume that the Zeeman energy is greater than the 
quadrupole energy, so that the nucleus possesses four nonequidistant spin energy 
levels E m with eigenfunctions % m ( m = ±3/2, ±1/2), which are eigenfunctions of the z 
component I z of the nuclear spin I. 

Let a radio frequency (RF) field 2H a cosQ a t, polarized along the y axis and 
resonant for a certain pair of energy levels M2 a = E m - E n , act on the spin. The 
interaction operator H a with the field in the interaction representation contains a time- 
independent term H a> eff = #yH a [<ml/ Y ln>P mn + <nl/ y lm> P nm ] and rapidly oscillating 

(at frequencies which are sums of Q a and (E k - Ei)//z) terms. The role of these last 
terms in the evolution of the spin states is negligible. To simplify the notation here 
and below we use the projection operators P mn representation , where P mn is a 4x4 
matrix with all matrix elements pki equal to zero except p mn =1. Projection operators 
are very convenient because of their very simple properties: PkiP m n = SimPkn and P mn %k 
= 5 n k%m- In addition, 1,2,3, and 4 will stand for the indices -3/2, -1/2, +1/2, +3/2, 
respectively. 

The interaction H a>eff leads to spin states evolution managed by propagator 

U(<p a ) = [P mm + P nn ] Cos((p a /2) + [P kk + P]J + [Pnm " Pmn] Sin((p a /2), (1) 



where (p a = 2(y H a t)kml/ ln>l, E m > E n and the indices k,l ^ m,n. If the RF field 
2HfCosn mn t is polarized along the x axis, the evolution operator is 

Uf((pf) = [Pmm + Pnn ] Cos((p f /2) + [P kk + P u ] + i [P nm + P ran ] Sin(<p f /2), (la) 
where (p f = 2(yH f t)kml/ x ln>l. 

In the currently accepted NMR model of a quantum computer two r<?<a/ 
exchange coupled spins R=l/2 and S=l/2 are considered as a basis for constructing 
quantum logic elements [1-4]. In the formalism of quantum mechanics the states of 
such a system and the operations on them are written in abstract four dimensional 
space, which is a direct product Er®Fs of the two-dimensional spaces of eigenstates 
of the real spins R and S. In our case, to clarify the information aspect of the proposed 
logical operations it is convenient to perform the inverse procedure: to represent the 
four-dimensional space T h corresponding to real spin 3/2, as a direct product r R ®r s 
of two abstract two-dimensional states spaces of virtual spins R and S. Then any 
operator P in the four-dimensional basis can be expressed as a linear combination of 
the direct products R ® S of operators given in the subspaces T R and Fs- The 
following isomorphic correspondence exists between the basis lm> of the space T t and 
the basis Im >®lm > of the direct product Tr®!^: 

I%i>oc |-l/2 > 1-1/2 > = II 1>, l% 3 > « 1+1/2 > I -1/2 > = I01>, 

\% 2 > I -1/2 > I +1/2 > = I10>, l% 4 > « I +1/2 > I +1/2 > = I00>, 

where I10> and so on are the notations adopted in information theory for states of two 
quantum bits (qubits). The energies corresponding to these states satisfy conditions: 
Ei>E 2 >E 3 >E 4 . 

The initial state for known quantum algorithms suggested for an abstract 
quantum computer is the state I00><*l%4>- From the standpoint of the quantum 
calculations the subsequent spin density matrix 



Pinit = const [ 1 + const P 44 ], 



(2) 



is the equivalent of the state I00>. Here 1 is the unit matrix in the space Tj. It does not 
change under unitary computational transformations and does not contribute to the 
observed signal when the result is read out. 

In a sample of macroscopic size a collection of such nuclei forms an ensemble, 
whose spin levels in equilibrium are populated according to the Boltzmann 
distribution 

p eq = Z exp( - (3 H ), Z l = Sp[exp( - (3 H )], p = 1/kT, (3) 

and at room temperature the relative populations difference of the spin levels is of the 
order of 10" 5 - 10" 6 or less. Therefore, to obtain the state Pi m t=constP 4 4 directly by 
cooling one requires very low temperatures, which, besides substantial technological 
difficulties, will slow down the speed of the entire computational cycle. 

To obtain the desired initial state for calculation we propose the procedure 
whose idea goes back to paper [3]. Let the required calculation consist of performing 
the transformation U CO mp of the state p m i t , while the spin system is in a state given by 
the density matrix (3), which in the high-temperature approximation is 

p eq = Z[l+ZX m P mm ], 1= £P mm (4) 

where X m = E m / kT and m = 1,2,3,4. We shall show that the sum of the results of 

three specially designed transformations of the state p eq is equivalent to a 
transformation of the state p init . Indeed, assume that we are required to perform a 
calculation which is a unitary transformation U com p of the spin I states. We define the 
unitary transformations 

Ui = U a (7U)U b (7U) = P13 + P 2 1+ P32+P44, 
U 2 = U b (7U)U a (7U) = - P12 - P23+ P31+P44, 

where U a (7t) is the propagator (1) for a pulsed RF field at frequency Q. d = (Ei - E2)/fi, 
which performs a rotation by the angle (p = %, and Ub(7t) is the propagator for a 
resonant RF pulse at frequency £lb = (E2 - E^/h. It can be verified that the average of 
the three transformations U comp , U com pUi and U comp U2, 



(l/3)[U com pP e q U comp U com pUi p e q U lU CO mp U com pU2 Peq U 2 U C omp ] 
=(l/3)U CO mp( Peq + Ul p e q U + l + U2 Peq U + 2 ) U + C omp =U C ompPinit U + comp 

is equivalent to the transformation U com p of the density matrix 

p imt = Z [al + (3P 44 ] , (5) 

where a = l+(l/3)[Xi +X 2 + X 3 ], and (3 = X 4 - (1/3) [Xi +X 2 + X 3 ]. 

Let us find the propagator corresponding to the rotation by a certain angle in 
the space T s under the condition that the space T R is invariant. This can be done by 
acting on the spin I with a two-frequency RF pulse containing the resonance 
frequencies Q.. d = (Ei - E2)/fi and £l c = (E3- E 4 )/fi. The propagator of such a 
transformation in the space Ti will be the product of propagators of the form (1) which 
perform rotation by the same angle (p a = (p c = (p at each transition 

U a+C (cp, q>) =1 Cos((p/2) + [P21 - P12 + P43 - P34] Sin(<p/2). 
It can be expressed as follows in terms of the operators of the spaces T R and Ts: 

U a+C (cp, 9) = (Rn + R 22 )®[ (Sn + S 22 ) Cos(<p/2) + (S21 - Si 2 )Sin((p/2)] = 

= exp{i((p/2)l R ®S Y }, (6) 
which proves the assertion made above (for virtual spins the indices 1 and 2 are used 
instead of -1/2 and +1/2, respectively; unit matrices are defined as 1 R = H R m m an d 

Is = £ S mm , with m = 1,2). Similarly, the propagator for a two-frequency RF pulse 

with carrier (filling) frequencies Q.^ = (Ei - E 3 )/# and Q. e = (E 2 - E 4 )/fi and angles (pa = 
(p e = (p will be equal to 

U d+e (<p,q>) = [P22 + P44 ] Cos(<p/2) + [P 42 - P 24 ] Sin((p/2) + 
+ [P 33 + P11 ] Cos((p/2) + [P31 - P13] Sin((p/2). 
It can be expressed in terms of spaces T R and T s spin operators as 

U d+e (<p,<p) = [ (Rn + R 22 ) Cos(<p/2) + (R 2] - R i 2 )Sin((p/2)] ® (Sn + S 22 ) = 



= exp{i((p/2)R x ®l s }. (7) 

which is the rotation by the angle (p in the space T R with Ts remaining unchanged. 
Generally speaking, the transitions at the frequencies Q. d and Q. e are forbidden in the 
initially adopted configuration of the magnetic and crystalline electric fields. Here we 
shall assume that gradient electric field asymmetry parameter is not equal zero and a 
small entanglement of the wave functions takes place. This allows «forbidden» 
transitions, and a larger amplitude of the RF field provides the required rotation angle. 
All results of the present letter can be obtained also with an arbitrary configuration of 
static magnetic and electric fields, but this would unjustifiably complicate the paper. 
Next, the transformation Uf ( (pf) at (pf = 71, defined as 

Uf(7C)= [P33+P44H i[P 2 i + Pi 2 ], 

performs the two-bit operation "controlled negation" CNOT - it performs the 
operation NOT on spin S if the spin R is in state ll> and leaves the spin S unchanged 
if the spin R is in state I0>. Indeed, it is easy to check that 

Uf(7C)l%i> =U t {7C) lll> =I10>, Uf(7C)l% 2 > =U t <7C) I10> =lll>, 

Uf(7C) l% 3 > = Uf(7C) I01> = I01>, Uf(7C) l% 4 > = Uf(7C) I00> = I00>. 

Hence one can see that the evolution operator Uf(7l) can be represented in the basis 

r R ®r s as 

Uf(7t) = 10x01 ®l s + llxll ® S x = R 22 ®ls+ R 11 ® S x , (8) 

which proves the assertion made above. 

To find the result of the computations it is necessary to read out the state of the 
final density matrix p ou t = U CO mp Pinit U + CO m P . NMR methods make it possible to 
measure all elements of the density matrix by the state tomography [4]. As an 
illustration, we shall examine a variant of reading out in the particular case of a 
diagonal density matrix 



pout=M +^(I m Pmm, m=l,2,3,4 



(9) 



in the important situation where the result of the computation is one of the states 
I Xm>, i-e-, when only one of the quantities |i m is nonzero. It is necessary to affect the 
real spin 3/2 with a pulsed two-frequency RF field, which generates a free-precession 
signal at the resonance frequencies Q. y2 and Q 34 by rotating the density matrix 

elements by the angles cp a =cp c =7l/2. Such a pulse corresponds to the evolution operator 

U 3 = U a (7u/2)U c (7t/2) = (1/V2)[ 1 +P 21 - P 12 + P 43 - P34] (10) 

Under the evolution operator (10) the density matrix (9) in the Schrodinger 
representation takes the form 

p(t)=(l/2){ ((X x +^ 2 )(Pn+ P22) + (U4+1I3XP33+ P 44 ) + 

+ ((1 1 - (I 2 )[ P 2 iexp(-itQ n ) + Pi 2 exp(itC2 12 )] + (11) 

+((j, 3 - M- 4)[ P43 exp(-itC2 34 ) + P34 exp(itft 34 )] }, 

where time is measured from the end of the computational cycle. In the state described 
by the density matrix (11), the quantum-mechanical averages of the transverse spin 
components become nonzero 

< I + (t)> = <I x +iI y > = <V3(P 43 +P 2 i) +2 P 32 > = Sp[p(t)(I x +iIy)] = 

= V3(|i 3 - [I 4) exp(-itQ 34 ) + V3((X 1 - (I 2) exp(-itQ 12 ), (12) 

The precession of the nuclear spin placed in the detection coil generates a periodic 
voltage at the two resonance frequencies Q. n and Q 34 with Fourier components V3((l 3 

- (i 4 ) and V3((ii - (i 2 ). We note that an identical pulse acting on the equilibrium 
density matrix (4) would produce a similar precession but with Fourier components 
equal to V3 Z(X 3 - X 4 ) and V3 Z(ki - X 2 ). A measurement of the sign of the ratios 

b 34 = (|i 3 - |i 4 )/ (A, 3 - A, 4 ) and b i2 = (|ii - (i 2 )/(^i - A, 2 ) 

of the corresponding Fourier components before and after the computation permits 
determining the final state of two virtual spins: 



if b34 < and bn = 0, then the result of the computation is I00>, 

if b 34 > and bn = 0, then the result of the computation is I01>, 

if D34 = and bn < 0, then the result of the computation is 1 10>, 

if b 34 = and bn > 0, then the result of the computation is 1 1 1>. 

Nuclear spin 3/2 is not rare. Nuclei with such spin appear in the most divers 
and easily accessible materials. Writing of two qubits on discrete levels of a single 
quantum particle obviates the need for an exchange interaction between the carriers of 
the qubits. In existing schemes for implementing quantum gates one needs to apply 
the special methods in order to suppress these exchange interactions during certain 
time intervals. This complicate the implementation of algorithms. The scheme 
described above can be applied to quantum systems of arbitrary physical nature. In 
principle, the same purpose one can reached by using nuclear spins of large magnitude 
(after choosing four suitable energy levels), ESR spectra with effective spin S* > 3/2, 
and optical energy levels. Only the expressions for the resonance frequencies and the 
matrix elements of the operators will be changed. Specifically, one can use a nuclear 
quadrupole resonance energy spectrum, split (or not split) by the interaction with a 
static magnetic field. The use of effective spins S* > 3/2 and therefore a large number 
of discrete energy levels gives, generally speaking, additional possibilities, which are 
not discussed in the present work. 
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